Short note on the Rabi model 
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The spectral density of the Rabi model is calculated exactly within a continued fraction approach. 
It is shown that the method yields a convergent solution. 
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A recent work by Braak [l| has renewed the interest in the old problem of coupling a photon field 
to a single spin 1/2 state, using the Rabi model [2]. The central statement of this work is that the 
eigenfunctions in Bargmann representation must be analytic functions in the entire complex plane. Based 
on this condition, a procedure is derived from the series expansion of the eigenstates which provides a 
recursive evaluation of the spectrum. It is claimed that the series expansion yields an exact solution of the 
Rabi model, which cannot be obtained from a direct continued fraction (CF) approach. In the following 
it is shown that this statement is not valid, and that the use of the extra condition of analyticity of the 
eigenfunction in Bargmann representation is not necessary. Therefore, the CF is directly applicable to 
the Rabi model and yields an exact solution for the spectral density, where the term "exact" means that 
the CF is convergent and the spectral density can be evaluated by a simple algorithm to any desired 
accuracy. 

The Rabi model is defined by the Hamiltonian 

Hji = cua^a + Ao- 3 + U(a 1i +a)cT 1 , (1) 

where (a) are creation (annihilation) operators of a photon, the Pauli matrices <Xj (j = 0,...,3) 
describe operations on the spin-1/2 state, and A is a symmetry breaking field for the spin. U is the 
coupling strength between the photons and the spin states. This Hamiltonian maps a product state 
\N) (g> |cr) with N photons to \N ± 1) <g> \a'}, where a, a 1 = , f- Thus eigenstates of H are superpositions 
of product states, which can be derived by a recursive approach for the coefficients of the superposition 
Q. In the following we apply the recursive projection method (RPM) of Ref. Q. This method is based 
on a random walk expansion in the underlying Hilbert space which visits each subspace at most once § . 
Consequently, there are no loops in the random walk, which leads directly to a CF. 

The RPM, as described in Ref. [4], can be directly applied to the resolvent (N; s\(z — H)~ 1 \N; s'} and 
yields immediately the spectral density with respect to \N; s): 

Pn{E) = \ Y Urn Im({N; s\(E - ie - H)- 1 ^; s)) . (2) 

s=t,l 

The simplest case is the resolvent (0;s\(z — -ff) _1 |0;s') of states without photons |0, s) and with spin 
projection s =f, 

({0;s\(z-H)- l \0;s'))=lim o ^ , (3) 

where the eigenvalues of H are the poles of the resolvent. Then the matrix elements g n , h n are subject 
to the following recurrence relations 

9k = z - u{n — fc) + A - U 2 {n -k + l)h k ^ ^ 

(5) 



- uj(n - fc) - A - U 2 (n - k + l)g k -i 
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for 1 < k < n with the initial values 

.90 = ; t- , h = , a • ( 6 ) 

z — urn — A z — uin + A 

The iteration of these recurrence relations gives a finite CF of the standard form Q 



9n = o 



Oil , a\\ a n \ 

= b o + TT- + ■ ■ ■ + tt 1 (7) 



with coefhcients 

oi = l, a k = (l-k)U 2 , o = 0, bi=z-A, h = z - (k - l)u> + (-l) k A (2 < k < n) . (8) 

For the corresponding coefhcients of h n we must only replace A by —A. Therefore, all considerations for 
g n apply to h n after replacing A — > — A. Due to b Q — the inverse of g n can be obtained directly by 
inverting g n in Eq. ([7]): 

a,/gn = b 1 + ^ + ... + ^. (9) 

This means that the poles of g n are obtained from the zeros of the related CF, and vice versa. 

By taking the limit g = lim n _^oo g-n the CF in Eq. (0 or ((9J yields the spectral density of the Rabi 
model for the complex energy z. This requires the limit n — >■ oo though, whose existence is a consequence 
of the Pringsheim Theorem (cf. Considering the tail t„ of the CF 

5 = °o + Tr + -- - + T^+*n ( 10 ) 
\0\ \o n 

with 

*n = ^ + -, (ID 
l&n+l 

we estimate the error when we truncate g at k — n (which gives g n ). The Pringsheim Theorem states 
that for complex a k , bk with \bk\ > \a.k \ + 1 (k > n) the truncated CF t n ,i 

+ a n+l| , a i| / . 1 . n / 10 , 

= ti 1 V TT (n + l<l) (12) 

\o n +i \t>i 

is convergent for I — >■ oo and converges to a t n with 

0<|t„|<l. (13) 

For the coefficients of the Rabi model in Eq. (jHJ the conditions for the coefficients are satisfied for a given 
2 and for sufficiently large n if U 2 < cj. Then (p~3|) means that the ratio t n /b n is 0(l/n), which implies 

This result enables us to approach the exact solution with a sequence of truncated CF's with a given 
error 0(l/n): 

, , a l\, , gn-jj , a n /b n + 0(l/n)\ 

9 = bo + jj- + ■ • • + TT + n ■ (15) 

\bi \b n -i |1 

The finite CF can also be expressed as the ratio 

= A, t _i + [a n /b n + 0(l/ra)]A n _ 2 
5 B n _i + [a„/6„ + 0(l/n)]B„_ 2 



(16) 
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FIG. 1: Spectral density po(E) for maximally n = 500 photons and A = 0.4, U = 0.99, e = 0.005 in units of the 
photon frequency ui. The green (red) curves are contributions from g„ {h n ). 



where A n _i (-EL-/) are polynomials in z of order n — I — 1 (n — I) which are generated by the recurrence 
relations (cf. 6]) 

A k = b k A k -i + a k A k _ 2 , B k = b k B k - X + a k B k _ 2 (k > 1) (17) 

with initial conditions A_i = 1, Aq = bo = 0, B-\ = 0, Bq = 1. 

Conclusions and remarks: Eq. (|16j) allows us to approach successively the spectral density of the Rabi 
Hamiltonian, where the accuracy is improved with increasing n. The error is estimated by the tail t n of 
(IT4l) . According to our definition this provides an exact solution of the Rabi spectrum, where the poles of 
g and h are the corresponding eigenvalues. The evaluation of pn{E) for N > can be performed within 
the same approach. 

The individual matrix elements g, h avoid level crossing due to parity conservation, since eigenstates 
of consecutive eigenvalues have different parity. However, the sum of g n and h n in the spectral density 
Po(E) provides level crossing by shifting their levels relative to each other, for instance, by changing A. 

The "rotating- wave" approximation of the Rabi Hamiltonian yields the Jaynes-Cummings Hamiltonian 

Hj C = wata + A<r 3 + U(a^a- + aa+) (18) 

with a± = (a i ±i<72)/2. This Hamiltonian does not belong to the class of models with recursive solutions 
because the recurrence relation of the RPM terminates already after the first step for each k: 

z- uj(n- k + 1) - A 
9k = [z - uj{n -k + l)] 2 - A 2 +cu - U 2 (n - k + 1) 



h k = 



1 

z — ui(n — k) 



(19) 



Here the evaluation of the eigenvalues requires onh 
agreement with the original work on this model 0, f 



the solution of a quadratic equation for each k, in 
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